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A simple analytic model of point-ion electrostatics has been previously proposed (R. Magri, S. -H. 
Wei, and A. Zunger, Phys. Rev. B 42, 11388 (1990)) in which the magnitude of the net charge qi on 
each atom in an ordered or random alloy depends linearly on the number N± of unlike neighbors in 
its first coordination shell. Point charges extracted from recent large supercell (256-432 atom) local 
density approximation (LDA) calculations of Cui-^Zn^ random alloys now enable an assessment of 
the physical validity and accuracy of the simple model. We find that this model accurately describes 
(i) the trends in qi vs. particularly for fee alloys, (ii) the magnitudes of total electrostatic 

energies in random alloys, (iii) the relationships between constant-occupation-averaged charges {qi) 
and Coulomb shifts (Vi) (i.e., the average over all sites occupied by either A or B atoms) in the 
random alloy, and (iv) the linear relation between the site charge qi and the constant-charge-averaged 
Coulomb shift Vi (i.e., the average over all sites with the same charge) for fee alloys. However, for 
bec alloys the fluctuations predicted by the model in the qi vs Vi relation exceed those found in the 
LDA supercell calculations. We find that (a) the fluctuations present in the model have a vanishing 
contribution to the electrostatic energy, (b) Generalizing the model to include a dependence of the 
charge on the atoms in the first three (two) shells in bec (fee) - rather than the first shell only 
- removes the fluctuations, in complete agreement with the LDA data. We also demonstrate an 
efficient way to extract charge transfer parameters of the generalized model from LDA calculations 
on small unit cells. 



PACS numbers: 61.66.Dk, 71.10.-fx, 61.50.Lt 



I. INTRODUCTION 

The structural stability of alloys and compounds is 
determined by the kinetic, electrostatic, and exchange- 
correlation contributions to the total energy. In first- 
principles calculations based on Hartree-Fock or on den- 
sity functional theory, the electrostatic portion of the 
total energy is characterized in terms of the electronic 
charge density p(r) and the nuclear charges 2^. For sys- 
tems with uniquely specified nuclear positions {Ri} and 
charges {zi}, the charge density is a well defined quantity 
as is the electrostatic (el) portion of the total energy: 
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structuEe^ffl'D'Elfl'LfLl and "supercell" models of random 
alloys|]'t3 there are well-defined {R^z^}, so the elec- 
trostatic energy was oMaiBftd-jpjm Eq. (|l|). However, 
in simpler approaches ,l 11 H 12 H ls lr 4 H 15 l one approximates the 
electrostatic energy by replacing the continuous charge 
density p(r), with fictitious point charges qi at each site i. 
For a system with N sites, the electrostatic or Madelung 
(M) energy is 
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where Rij is the distance between sites i and j. The 
Madelung energy may also be written 
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where Vi is the Coulomb shift at site i due to all charges 
other than q^. 
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The point charges are obtained by partitioning p(r) into 
"domains" (spheres, polyhedra, etc.) and integrating the 
total charge in each domain. However, because there is 
not a unique way to partition a three-dimensional space, 
the point charges are not uniquely defined. 

For periodic systems (e.g., ordered structures with 
a primitive, .( jell . pij |-jandom structures defined by 
swperce/^BilJ'OllaO'ta) where all sites i are defined as 
distinct entities (not as averages) and qi and Rj are spec- 
ified, Em can be readily computed from Eq. (|J) using, 
for example, the Ewald method. In most statistical ap- 
proaches to alloy | s-.(e.g., the coherent potential approxi- 
mation, or CPA)E3 however, one attempts a description 
of a random alloy without a specification of all distinct 
sites i but rather some averages over i. In such ap- 
proaches one calculates the Madelung energy of the ran- 
dom alloy by determining the configurationally averaged 
correlation between charges (qiqj), and using 
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Until 1990, all CPA-based models for alloy energies have 
assumed uncorrelated charges 



(QiQj) = 



(0) 



which leads to a vanishing electrostatic energy for the 
random alloy 
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on account of electroneutrality. This approximation [Eq. 
(0)] was based on the expectation that a random (i.e. 
uncorrelated) distribution of atoms on sites would lead 
to an equally random distribution of charges, i.e., the 
charge on an atom in a given alloy is a property of the 
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atom, irrespective of its environment. Eq 
assumed in many CPA-based calculations!! 
involving the total energy of random alloys. Magri et 
aZ.EJ subsequently criticized this approach as being phys- 
ically implausible, since the assumption of uncorrelated 
charges [Eq. (||)] means that an A atom surrounded lo- 
cally by only A atoms will have the same charge as an 
A atom surrounded by B atoms; chemical intuition sug- 
gests, however, that the charge on a site will depend on 
the identity of atoms in its environment because charge 
transfer is present only between dissimilar sites. 

Magri et alx3 noted that in a random alloy, even 
though the occupation of site i is independent of the occu- 
pation of other sites by definition, the charges on a site do 
depend on the occupations of other sites. These authors 
therefore proposed a simple model to describe the magni- 
tude of point charges in disordered (and ordered) alloys: 
The magnitude of the charge on a site is linearly pro- 
portional to the number iv| of unlike nearest neighbors 
surrounding that site. With this charge model, Magri et 
al. went on to demonstrate that even for the case of a 
random alloy with completely uncorrelated atomic occu- 
pations, charge correlations exist in the alloy and these 
correlations lead to a non-zero Madelung energy. 

Subsequent to the proposal of Magri et a/. jLil the charge 
model has been-used in many contexts: 

(i) Lu et alS3 showed that LDA calculations on or- 
dered compounds produced charge densities which, when 
integrated inside muffin-tin spheres, gave point charges 
which reproduced the behavior of the model. They also 
examined! the effect of the ensuing electrostatic energy 
of the random alloy on-the sign of the ordering energy. 

(ii) Abrikosov et alt5 and Johnson and Pinskiliij de- 
rived corrections to the CPA total energy which intro- 
duced charge correlations in random alloys. These cor- 
rections were shown to be consistent with the charge 
model of Magri et al. Several authors subsequently 
used these corrections in total-energy CPA calculations 
to determine lattice constants and formation energies 
of random metallic alloys, finding significant effects 
due to charge correlations: Johnson and Pinskit^l es- 
timated the total energy contribution due to charge 
correlations to be —1.25, —5.3, and —7.7 mRy/atom 
for CU0.5Z110.5, CU0.5AU0.5 and Nio.sAlo.s alloys, respec- 
tively. (Typical values of alloy formation energies are 



^10-20 mRy/atom.) Korzhavyi et al. founded that 
the energetic contribution due to charge correlations for 
Alo.sLio.s is -16.0 mRy/atom which results in a change 
of sign in the formation enetgy of Al-Li alloys. 

(iii) Borici and Monnierca used the charge model to 
study the segregation behavior of a semi-infinite ran- 
dom Madelung lattice. For semi-infinite surface geome- 
tries, these authors found that charge correlations lead 
to monotonic surface segregation profiles and a segrega- 
tion of the minority species to the surface. On the other 
hand, for thin-film geometries charge correlations lead to 
oscillatory surface segregation profile, and an enrichment 
of the majority species on thej-sprface. 

(iv) Wolverton and Zungeili-3 determined the ground 
state long-range order and the high-temperature short- 
range order of fee-, bec-, and sc-based alloySpdue to elec- 
trostatic effects. These authors also showedElo how the 
charge model could be analytically mapped onto a cluster 
expansion, which allowed for the efficient and accurate 
determination of energies of any ordered or disordered 
configuration without the use of Ewald methods. 

(v) Ruban et al£3 compared the energies of charge- 
correlated CPA calculations with ordered compound 
LDA calculations to determine the optimum prefactor 
for the electrostatic energy for Cu-Au and Ni-Pt alloys. 
The energetic contribution due to charge correlations was 
again found to be significant: For instance, for random 
CU75AU25 alloys, electrostatic contributions to the total 
energy were found to lower the mixing energy by a fac- 
tor of ~3-6 relative to CPA calculations with a complete 
neglect of charge transfer effects. 

The charge model ansatz of Magri et al. was thus far 
tested by comparing its charges {qi} with those found 
in small-unit-cell (<16 atom) LDA calculations, and also 
for only one lattice type - fee. Recently, jhucIi larger 
LDA supercell calculations became availableMI for fee and 
bec-based alloys. These calculations combine a locally 
self-consistent muffin-tin scheme with a massively paral- 
lel computer enabling LDA calculations o:a-2j56.- and 432- 
atom-aunercells for random Cu-Zn alloysOEa Faulkner 
et a^liZEj have used the charge density from these large 
LDA supercells to examine the behavior of point charges 
{qi} in random Cu-Zn alloys, finding interesting relations 
between charges and certain potentials. Here we deter- 
mine to what extent the simple charge model is able to 
describe the electrostatic properties of complicated large 
scale (256-432 atom) LDA based calculations. We find 
that the model works very well for fee lattices, but that 
in bec lattices, where the first few coordination shells are 
near to one another, the charge on a site is correlated with 
the occupations on a few neighboring shells, not just one. 
The effects of such corrections to the total electrostatic 
energy (Em)r are small, however. 



II. THE SIMPLE CHARGE MODEL 

Consider an A\^ X B X alloy with N sites and a nearest 
neighbor coordination number Z . The model of Magri et 
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is based on the assumption that the excess charge on 
a site depends only on the identity of its first neighbors. 
If an A atom on a central site is surrounded purely by Z 
atoms of type A, the charge is taken to be zero. If it is 
surrounded by Z atoms of type B, the charge is maximal, 
2ZX. For intermediate occupations of the first coordina- 
tion shell, we assume a linear interpolation between these 
two limits. Formally, we then write this charge as: 



q% = A ^[Si - Si+k], 



(8) 
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where the pseudospin Si is -1 (+1) if an A(B) atom is 
located at site i. (The set of variables Si for all sites i de- 
fines the configuration a.) Si+k indicates the occupation 
of the Z lattice sites which are nearest neighbors to i, 
and hence the summation in Eq. (|^) indicates the num- 
ber of unlike nearest neighbors surrounding the site i. X 
is a constant which indicates the magnitude of the charge 
transfer and is an undetermined parameter of the model. 
Thus, the charge model will give trends in the behavior 
of physical properties, but will not give numerical values 
of properties without some input value of A. 

Several questions may be asked concerning the param- 
eter A: 1) Should A be explicitly composition-dependent? 
2) Should A be explicitly volume-dependent? Since the 
equilibrium volume is a function of composition in size- 
mismatched alloys, an explicit volume dependence of A 
would lead to an implicit composition-dependence. It 
is important to physically distinguish between these two 
dependences. 3) Should the values of A be extracted from 
large unit cell or small unit-cell alloys i.e., does A con- 
tain mostly short-range or long-range information ? .. VaL 
ues of A have been estimated by LDA calculations ,E3£3'E£l 
ranging from small unit cell ordered compounds (~8-16 
atoms), up to large LDA simulations of random alloys 
(^200-400 atoms). For computational simplicity, one 
should know whether it is equally valid to extract val- 
ues of A from ordered or random alloys, and whether one 
can even use smaller cells (~2-4 atoms) than have been 
currently used. 

We next examine the physical consequences of charges 
which obey Eq. dsl). We then compare these conse- 
quences with results of LDA supercell calculations in or- 
der to assess the physical validity of the model. With re- 
gard to the questions raised above, we demonstrate that 
the simple charge model represents well the charge trans- 
fer of different unrelaxed configurations at a common vol- 
ume. If more than one volume is considered (e.g., for a 
lattice- mismatched alloy at more than one composition) , 
the parameter of the model A would presumably need to 
be explicitly volume-dependent (implicitly composition- 
dependent). Also, we find that values of A extracted from 
2-4 atom LDA calculations agree favorably with those 
extracted from much larger 200-400 atom LDA calcula- 
tions, thereby resulting in a drastic computational sim- 
plification. 



III. PHYSICAL CONSEQUENCES OF THE 
CHARGE MODEL 

A. Average Charges 

The average charge on all sites, (q) is defined as: 

i 

Combining this with Eq. (||) gives (q) = 0, as guaranteed 
by global charge neutrality. However, what is more inter- 
esting is the constant-occupation- average (q) a (or (g)s), 
i.e., the average charge of all sites occupied by A (B) 
atoms. This constant-occupation-average is a function 
of the configuration a and composition x, and we can 
analytically derive this quantity for any arbitrary config- 
uration. The definition of {q) a is 
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where Na is the number of A atoms in a and is the 
Flinn operator such that Tf = 1 if site i is occupied by 
an A atom, and Tf — otherwise. The Flinn operator 
is given by Tf = (l- S l )/2. Thus, 
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where II is the nearest neighbor (NN) pair correlation 
function, i.e., the lattice average of = SiSj for i and 
j NN. A similar analysis gives 
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In addition, the difference A in constant-occupation- 
averaged charges is given by 



A = (q) B - (q)a = 2ZX 



i -n 



(13) 



where (S) = 2x - 1. Equations Q, @, and |[|) are 
quite general and apply to any configuration (ordered, 
random, or partially ordered). These expressions may be 
evaluated in various classes of configurations which are 
interesting: 

Random Alloys: We define a random alloy as one in 
which the occupation variables Si are uncorrelated, i.e. 
IT j = (Si)(Sj) = (2x — l) 2 . As we will see below, this 
does not imply that the charges qi are uncorrelated. In 
a random alloy, we have 



-2ZXx ; (q) B = 2ZX(1 - x) ; A = 2ZX (14) 



Interestingly, A for random alloys is independent of alloy 
composition x. 
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Short-range Ordered Alloys: Short-range order (SRO) 
measures the extent to which there are atom-atom cor- 
relations in a disordered alloy. The relation between the 
nearest neighbor Warren-Cowley SRO parameter (a) and 

n is 



n-(s) 2 

(5)2 
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This expression combined with Eqs. ([Tl]) and (|lj) and 
( |l3| ) gives the constant-occupation-averaged charges for 
alloys possessing some degree of SRO: 

(q) A = -2ZXx(l-a) ; (q) B = 2Z\(l - x)(l - a) ; 

A = 2ZX(l-a) (16) 

So, the difference in charges A should increase in an or- 
dering type alloy (a < 0) relative to the random values, 
but should decrease in a clustering type alloy (a > 0). 

Long-range Ordered Alloys: Long-range order (LRO) 
gives an indication for the relative population of A or B 
atoms on a given sublattice. The extent of LRO in an 
alloy may be described by (one or more) LRO parameters 
n. For example, for an alloy at x=l/2 with a single LRO 
parameter < r\ < 1 (and no correlations between atoms 
on the same sublattice), H(rj) = ?7 2 II(1), so for any state 
of LRO at x—1/2: 



(q) A (ri)=-ZX[l- V y(l)] 
(«> B (t/)=ZA[l-rr|n(l)] ; 
Afa) = 2ZX[1 - n 2 n(l)} 



(17) 



For example, in CsCl (B2) type ordering 11(1) = —1, 
thus as the degree of LRO increases, A increases due to 
the increased number of unlike nearest neighbors. 



B. Constant-Occupation- Averaged Coulomb Shifts, 

(V)a and (V)b 

The Coulomb shift Vi [Eq. (^)] averaged over all sites 
(V) is zero (just as (q) is zero) due to global neutrality. 
A more interesting quantity is the constant-occupation- 
average of the Coulomb shifts on all sites occupied by 
A- atoms in a random alloy: 



(18) 



where the latter brackets denote a configurational aver- 
age. This expression can be evaluated to give: 



X 



2(1 -x) 



jjti k 



X 



3^1 k J 

2xZX 



j+k 



Ri 



(19) 



where R\ is the nearest-neighbor distance and in the sec- 
ond equality of Eq. (|l9|), we have used the orthonormal 
properties of the products of Si £3 Similarly for B-atoms, 
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C. Relation between Constant-Occupation- Averaged 
Charge and Coulomb Shift 

From Eqs. ( |l4|) and ([l9|), we have the following re- 
lation between constant-occupation-averaged charge and 
Coulomb shift, as predicted by the charge model: 
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D. Charge-Charge Correlation Functions 

The charge-charge correlation function between sites i 
and j is given by: 



(QiQj) = 9i+mQj+ 
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The sums over k (fc') are over the nearest neighbors of 
site i (j). Eq. ( |22[ ) is generally valid for any configuration 
and any composition. For random alloys (i.e., alloys with 
uncorrelated site occupations), the charge correlaticais for 
the mth shell (<Zo9m) have been previously derivedH and 
are given by: 



(<?o<?o) = 4x(l 
(<?o<?i) = 4x(l 
{qmm) = 4cc(l 



x)X 2 {Zl + Z 1 ) 
x)X 2 (-2Z 1 + K 1 ) 
x)X 2 {K m ) ; m > 1 



(23) 



In these expressions, Z m is the coordination of the mth 
shell (i.e., Z\ = Z), and K m is the number of nearest- 
neighbor atoms-shared by sites i and i + m. As found 
by Magri et alE3, Eqs. (|2^) demonstrate that even when 
the occupations of sites are uncorrelated, the charges on 
these sites, obeying Eq. (||), are correlated. 



E. Electrostatic Energies of Random Alloys 

Using the charge-charge correlations in Eq. (^3|), one 
can obtain the electrostatic energies of random alloys 
which are a consequence of the charge model. These 
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energies of random alloys have hepn derived previously 
for fee, bee, and sc-based alloys Jillo The energies of fee 
and bec-based random alloys are given by 

(EmYr/Eq = -4jc(1 - ar)0.7395181... 
(E M } h R C /E = -4.t(1 - a;)0.3457752... (24) 

where E = (16A) 2 /2i?i. 

F. The q — V Relation between Charges and 
Coulomb Shifts 

For a completely random alloy, we may analytically 
derive (Appendix A) a relation between the charges qi 
and Coulomb shifts Vi from the charge model: In the 
model of Eq. (||) , the magnitude of the charge on a site 

% surrounded by unlike neighbors in the hrst shell 

does not depend on the spatial configuration of the 

and is | qi [N^ } ] | = 2AA, (1) . The Coulomb shift Vi[N^} 

on a site surrounded by N* 1 ' unlike nearest neighbors on 
the other hand, does depend in the model of Eq. (||) 
on the spatial configuration of the unlike atoms 

around i and also depends on more distant neighbors. 
If we average over all sites having JV- unlike neighbors, 
we find analytically (Appendix A) the linear relation be- 
tween charge qi[N^] and the constant-charge-averaged 
Coulomb shift Vi[N^]: (i.e., an average over all sites 
with the same charge, and hence with the same N' t ■ ), 

qdN^^V^}, (25) 

where 7 (in Ry -1 ) is the slope of this linear relation: 

7fcc (a; = 1/2) = 0.132i?i 

7bcc(z = 1/2) = 0.16372! (26) 

Note that Vi is a constant-charge-average (still leav- 
ing the dependence), in contrast to the constant- 
occupation-averaged (V)a- To evaluate the fluctuations 

in Vi(Ni) about Vi, we perform large-unit-cell sim- 
ulations for a single, randomly selected configuration. 
Equal numbers of atomic types A and B (x=l/2) are 
distributed at random over the 256 fee sites and 432 bec 
sites of the simulation cell. Point charges {qi} are then as- 
signed by the model of Eq. (|J) . Using the Ewald method, 
we then calculate the Coulomb shifts Vi [Eq. (^)] for each 
site in the cell. This gives a qi vs. Vi relation for the 
charge model. 

IV. COMPARISON: SIMPLE MODEL VS. LDA 
SIMULATIONS 

Recently,0lil large scale LDA supercell calculations 
(256- and 432-atom) have been performed for Cu-Zn al- 
loys with Cu and Zn atoms placed randomly on the fee or 



bec lattice sites of the supercell. These calculations uti- 
lize a multiple scattering framework, and are locally self- 
consistent: The charge density associated with each atom 
is constructed by considering only the electronic multi- 
ple scattering processes in a finite spatial region (several 
neighboring shells) centered at that atom. These LDA 
calculations also use the muffin-tin approximation: The 
charge within the Wigner-Seitz cell surrounding each site 
i (volume tti) is made of a spherically symmetric portion 
p l (r) = p l MT {r) inside each muffin-tin sphere (r < R l MT ) 
and is equal to a constant po in the interstitial region 
between spheres. Point charges are then extracted from 
the muffin-tin charge density by performing the following 
integral: 

f fl MT 4^ 

Qi = 4tt y p l MT (r)r 2 dr + p [^i - -y(^MT) 3 ] - z i 

(27) 

The calculations are carried out for unrelaxed geome- 
tries, thus each of the 256 or 432 atoms has an equiv- 
alent Wigner-Seitz cell. Cu and Zn have a very small 
electronegativity difference, so Madelung energy is quite 
small in Cu-Zn alloys, and could be more susceptible 
to any errors in the calculation. Cu-Zn is therefore a 
critical test of any charge model, as the electrostatic ef- 
fects in this system are quite subtle. An alloy system 
with more robust charge transfer (larger electronegativ- 
ity difference) could therefore be of interest in comparing 
magnitudes of electrostatic energies, charges, and other 
properties. .— . 

Differences of ~ 2% are citedliZI between LDA calcula- 
tions using bec cells of 256 and 432 sites, and hence indi- 
cate a typical error due to using a finite-sized supercell. 
An additional consideration is that one, single configura- 
tion is considered rather than a configurational average. 
Thus, in comparing properties of the charge model with 
those of LDA supercell calculations, only disparities of 
more than a few percent should be considered meaning- 
ful. 



A. Dependence of Charges on the Nearest-neighbor 
Environment 

The prediction of the model of Eq. (@) for the depen- 
dence of charge on the number of unlike-nearest neigh- 
bors is clear: It is a linear relation. The charges pre- 
dicted by this simple model (open circles) are compared 
with those obtained in extensive LDA supercell calcula- 
tions (crosses)Ea in Fig. [|. A least-squares fit to the LDA 
supercell data gives values of the parameter A: 

A fcc = 0.00819, A bcc = 0.01176. (28) 

The LDA calculations demonstrate that (i) the linear pre- 
diction of the model is accurate for fee alloiZB—hut (ii) for 
bec alloys, as recognized by Faulkner et aZ.Ji3Ea there are 
fluctuations about the straight line. Recall that in an fee 
structure the distances from the origin to the n th shell 
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FIG. 1. Charge versus number of unlike nearest neighbors. 
Shown are the predictions of the charge model of Eq. (M) and 
the results of LDA supercell calculations of Ref. Values of 
A [Eq. (p^) ] were obtained by a least-squares fit to the LDA 
supercell data. 



supercell calculations are compared in Table Q. The sim- 
ple charge model of Eq. (||) reproduces R e g/R\ of the 
LDA supercell calculations (0.97 and 1.02 for fee and bec, 
respectively) to within a few percent. 

In Eq. (M), it is shown that the simple charge model 
predicts that the difference A in constant-occupation- 
averaged charges is independent of composition. The 
LDA supercell results (Fig. 3 of Ref. [l7]) also show that 
A depends very little on composition, in agreement with 
the model prediction. 

The difference between constant-occupation-averaged 
charges A (77) for the charge model was also derived as a 
function of LRO parameter in Eq. (p"7|). The charge 
model prediction is that A (77) should increase for the 
LRO B2 alloy relative to the random alloy bv-,a fac- 
tor A S2 (1)/A(0) = 2. The LDA calculations^ show 
that the introduction of LRO does increase A from a 
value of 0.20066 for the random alloy to 0.25178 for the 
B2 ordered alloy, giving A B2 (1)/A(0) = 1.25, somewhat 
smaller than but qualitatively consistent with the pre- 
diction of the charge model. In considering the dispar- 
ity between the magnitude of A B2 (1)/A(0) of the simple 
charge model and that of LDA, one should remember that 
in the latter, point charges are defined by a non-unique 
partitioning of space. 

The influence of SRO on A was derived for the charge 
model in Eq. ( |T6| ) where it was shown that ordering type 
SRO (as found in Cu-Zn alloys) should increase A relative 
to the random alloy. TheJxitroduction of SRO in the 
LDA supercell calculationsEZl increases A from 0.20066 
for a random simulation to 0.20554 in a simulation with 
some degree of SRO. This increase is again consistent 
with the predictions of the charge model. 



are R%, 1.414Ri, and 1.732i?i for n=l,2,3, while for bec 
these distances are 1.155i?i, and 1.633i?i. Thus, the 
bec structure shows a weaker distinction between first 
and second neighbors. We consider below (Section |v|) 
possible generalizations of the charge model for bee al- 
loys which account for these fluctuations by extending the 
linear relation (^) to more than one coordination shell. 

B. Relation between Constant-Occupation- Averaged 
Charge and Coulomb Shift 

Both the charge model of Eq. (||) and the LDA su- 
percell calculations result in a relationship between the 
constant-occupation-averaged charges {q)A and (q)s [Eq. 
(|l4|)l and the constant-occupation-averaged Coulomb 
shifts (V) A and (V) B [Eq. @] of the form: 

(V}a,b = (29) 

According to Eq. (^l[) , the charge model of Eq. (j|) gives 
R c ff = i?i (where R\ is the nearest neighbor distance). 
Values of R e g/Ri for the charge model and for the LDA 



C. The q — V Relation between Charges and 
Coulomb Shifts 

The large supercell LDA calculations find a linear re- 
lationship between the charges qi on individual sites and 
the Coulomb shifts Vi on those sites. The simple charge 
model of Eq. (||) predicts [Eq. |25| a linear behavior 
between charge and constant-charge-averaged Coulomb 
shift in agreement with the LDA supercell calculations. 
The slope 7 of this linear relation is compared with the 
slopes from the LDA supercell calculations in Table Q. 
(Note that both charge and Coulomb shift are propor- 
tional to the parameter A of the model, and thus, the 
slope 7 is independent of A.) The relative slope j/Ri of 
the model (-0.132) is within a few percent of the LDA 
supercell results for fee alloys (-0.125 ± 0.002), while for 
bee alloys, the slope of the model (-0.163) is too large in 
magnitude relative to the LDA result (-0.115 ± 0.001). 

The relationship between qi and the distinct Coulomb 
shift Vi (not the constant-charge-averaged Vi) as ob- 
tained in the simple model of Eq. (||) is shown in Fig. 
||, where it is contrasted with the results of the LDA cal- 
culations of Ref. [l^. The fluctuations in Coulomb shift 
about the average linear behavior of qi and Vi are quite 
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TABLE I. Comparison of physical properties of random alloys which are a consequence of (a) the charge model of Eq. (J8|), 
(b) the generalized charge model of Eq. (^) , and (c) those obtained from LDA supercell simulations of Refs. [l7| |l8| In cases 
(a) and (b), we assign charges to sites according to a given model [Eq. (Q) or (pS4|)] and then calculate the Coulomb shift Vi 
[Eq. mft] by applying the Ewald method to the assigned charges. 



Refi/Rl 



-r/Ri py" 



{E m )r [mRy/atom] 



fee (as=l/2) 



Model - Analytic 
Model - 256 atoms 
Model - 16 atom SQS 
LDA supercell - 256 atoms 
Gener. model - 256 atoms 
Gener. model - 16 atom SQS 



1.00 
1.06 
0.99 
0.97 
1.02 
1.03 



-0.132 
-0.130 
-0.134 
-0.123,-0.127 
-0.120 
-0.118 



-2.60 
-2.52 
-2.60 
-2.61 
-2.55 
-2.52 



bee (x=l/2) 



Model - Analytic 
Model - 432 atoms 
LDA supercell - 432 atoms 
Gener. Model - 432 atoms 



1.00 
0.98 
1.02 
1.19 



-0.163 
-0.155 
-0.114,-0.116 
-0.119 



-2.57 
-2.64 
-2.67 
-2.34 



small in the fee random alloy, but are substantial in the 
bee alloy. We have next determined the effect of these 
fluctuations on the electrostatic energy (Em)r of the ran- 
dom alloy: If the linear relation Eq. between charge 
and Coulomb shift (neglecting fluctuations) is used in 
Eq. (^) to compute the random alloy energy, we recover 
precisely the same energy including fluctuations derived 
in Eq. (|24|). Thus, although the fluctuations in Vi are 
graphically impressive (Fig. |2j) , the energetic consequence 
of these fluctuations is strictly zero, simply indicating 
that the fluctuations in Coulomb shift are symmetrical 
about the average linear behavior. 



for the more distant 3rd, 4th, and 5th shells are extremely 
sensitive to the size of the simulation cell: For a single 
fee 256-atom simulation, one can even find 3rd and 4th 
neighbor correlations which have opposite sign relative 
to the exact analytic values. Even for very large (16384- 
atom) fee simulations configurationally averaged over 20 
configurations, the 3rd and 4th neighbor correlations may 
differ from the analytic values by ~ 10%. Thus, in order 
to compare the LDA charge-charge correlations for ran- 
dom alloys with the analytic results of the simple charge 
model, the size of the LDA supercell calculations would 
have to be significantly increased. 



D. Charge-Charge Correlations Functions 

The simple charge model predicts specific values for the 
charge-charge correlations of random alloys given in Eq. 
(p3|). The quantitative results of the LDA_.supercell for 
these correlations are a bit unclear: In Ref£3, the authors 
note that for 256-atom LDA supercell calculations, the 
nearest-neighbor correlations are sizeable, but they also 
note that the values beyond the nearest-neighbor shell are 
smaller than the predictions of the model, although these 
values are not too well known due to the relatively small 
size of the simulation cell. When larger LDA supercell 
simulations become available, a comparison of charge- 
charge correlations from LDA supercell with the predic- 
tions of Eq. (^) (and those of the generalized charge 
model described below) would be of interest. The ana- 
lytic values of the charge-charge correlation functions of 
Eq. d23| ) are plotted in Fig. |^. We have also compared 
these analytic values with those obtained from our large- 
unit-cell simulations of the charge model (not shown). 
Although the correlations for the nearest neighbor shell 
are robust with respect to unit cell size, the correlations 



E. Coulomb Energy of Random Alloys 

The Madelung energy of the simple charge model for 
a random alloy is given in Eq. ( p4| ) in terms of the pa- 
rameter A and the nearest-neighbor distance R\. If we 
use the numerical values for A given in Eq. ( psj ) and 
the nearest neighbor distances used in the LDA supercell 
calculations for Cu-Zn alloys, 

R[ cc = 4.879a.u.; R\ cc = 4.763a.u., (30) 

we obtain the electrostatic energies of the simple charge 
model for Cu-Zn: 

(Em)r({cc] x = 0.5) = -2.60 mRy/atom, 
(E M ) R (icc; x = 0.7) = -2.18 mRy/atom, 
(E M ) R (bcc; x = 0.5) = -2.57 mRy/atom. (31) 

These-siames are compared with the Cu-Zn LDA supercell 
valued in Table |: 

(Em)r(^cc;x = 0.5) = -2.61 mRy/atom, 
(Em) fl(fcc; x = 0.7) = -2.20 mRy/atom, 
(-EW)fl(bcc; x = 0.5) = -2.67 mRy/atom. (32) 
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FIG. 2. Charges versus Coulomb Shifts as predicted by the charge model of Eq. (jS) using the values of A and _Ri given in 
Eqs. (pg|) and dio|). Also shown are the results of LDA supercell calculations of Ref. [l8[ 



For all cases, the prediction of the simple model is ex- 
tremely accurate: the model energies fall within 0.1 mRy 
of the LDA supercell calculations. Although the model 
of Eq. (|J) was shown to have significant fluctuations in 
the q — V relations (Fig. ^|), these fluctuations have a 
vanishing contribution to the Coulomb energy, and thus 
the model produces accurate energetics. 

F. Approximating Large Random Supercells by 
Small-Cell "Special Quasi-Random" Structures 

Our foregoing discussions were based on (either LDA or 
Ewald) simulations of rather large supercells (e.g., 256- 
432 atom). We next examine the extent to which spe- 
cially selected small cells can mimic larger, non-specialhi 
selected cells. Special quasi-random structures (SQS)lI 
are small unit-cell structures which are constructed in 
such a way so that structural (not charge-charge) corre- 
lation functions IIsqs match as closely as possible those 
of the random alloy (IIsqs ~ IlRandom) for several neigh- 
boring shells. In this way the SQS is a small-unit-cell 
ordered structure which mimics the random alloy. It is 
interesting to see how the charge model calculation of a 
small-unit-cell SQS compares with the large scale 256- 



atom simulations described above. 

We have performed an Ewald calculation for a 16- 
atom, fcc-based SQS structure (denoted SQS-16) with 
point charges taken from the charge model of Eq. (^). 
Structural information for SQS-16 is given in Appendix 
B. The resulting R c g/Ri, j/Ri, (Em)r, and the q — V 
relation for the SQS-16 are collected in Table |, where 
they are compared with analogous calculations using a 
(randomly selected) 256-atom cell. We see that the 16- 
atom SQS calculation matches the 256-atom simulation 
for all properties to within a few percent. Also, the elec- 
trostatic energy of the SQS-16 (E m )sqs/E = -0.740 
compares much more favorably with the exact value of 
-0.7395 than does the energy of the 256-atom simula- 
tion {E M ) 2 ™~ atom /E Q = -0.716. Thus, for the case of 
electrostatic energies, the 16-atom SQS provides a more 
accurate depiction of the random alloy than does a single, 
randomly-selected 256-atom configuration. It would be 
interesting to compare the LDA energies of the SQS-16 
with those of larger, but randomly selected supercells £3 
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FIG. 3. Charge-charge correlation functions for random x=l/2 alloys. Shown are the correlations of the simple charge model 
and the generalized charge model for fee and bee lattices. For the simple charge model, the analytic results of Eq. ( |23| ) are 
plotted while the correlations of the fee (bee) generalized charge model are obtained from a 16,384 (8,192)-atom simulation, 
configurationally averaged over 20 configurations. Correlation functions are shown normalized to Af 



V. GENERALIZATIONS OF THE MODEL 

A. Summary of Successes and Failures of the Simple 
Model 

We have thus far ascertained the physical predictions 
of the simple, nearest-neighbor charge model of Eq. (||), 
and compared them with the results of large LDA su- 
percell calculations of Ref. [l7|, [l8|. In many cases, the 
simple charge model accurately predicts the electrostatic 
properties of LDA: 

(i) The behavior of charge versus nearest neighbor en- 
vironment is reproduced well by LDA calculations of fec- 
based alloys (Fig. 0a). 

(ii) The proportionality R e g/Ri between constant- 
occupation-averaged charge {q) a.b and Coulomb shift 
(V)a,b of the model is the same as that of LDA to within 
a few percent (Table ||) . 

(iii) For fee alloys, the linear q — V relation of LDA 
is well reproduced (including fluctuations) by the model 
(Fig. ||). The model value for the slope 7/-R1 of the 
q — V relation in fee alloys is within a few percent of 
LDA (Table |). 

(iv) The Coulomb energies of the model are extremely 
accurate with respect to the LDA values (to within 0.1 
mRy/atom). 

(v) The slope of the q — V relation, 7, in Fig. || are the 
same for A and B. The LDA supercell calculations also 
show similar slopes (7/-R1) for charges on Cu (—0.123) 
or Zn (—0.127) atoms in the fee x—\/2 alloy, or for Cu 
(-0.114) or Zn (-0.116) atoms in the bec x=l/2 alloy. 

(vi) The slope of the charge versus number of unlike 
nearest neighbors (Fig. |l|) are negatives of one another. 
LDA supercell calculations (for fee alloys) support this 
(Fig- J). 

(vii) In the impurity limit, the model predicts that the 
charge on A embedded in pure B is equal (in magnitude) 



to that of B embedded in pure A, 

\{q) A (x - 1)| = \{q) B (x -» 0)| = 2ZX. (33) 

The LDA supercell calculations also show this behavior 
(see Fig. [|), for an atom surrounded completely by unlike 
neighbors. Note that neither the simple model nor the 
LDA supercell simulations include the effects of atomic 
relaxations, which could likely eliminate the degeneracy 
of Eq. (H|). [To describe relaxed configurations, it is 
anticipated that more parameters (e.g., bond lengths) 
would need to be introduced into the model.] 

(viii) A is composition-independent in the charge 
model; values of A (Table I of Ref.tLa) extracted from 
the LDA supercell calculations also demonstrate that A 
is not sensitive to concentration. We reiterate that the 
charge model describes only unrelaxed configurations at 
a fixed volume. For lattice-mismatched systems, alloys 
of different composition will have different volumes, and 
the charge transfer will depend on this volume. To model 
this effect, A should be explicitly volume-dependent. This 
explicit volume-dependence would lead to an implicit de- 
pendence of A on composition. (Presumably, this im- 
plicit composition-dependence is not seen in the LDA 
supercell data of RefJia due to the fact that the system 
studied, Cu-Zn, has a relatively small lattice-mismatch.) 
However, this should not be confused with an explicit 
composition-dependence of A. 

Although there are many cases of agreement between 
the predictions of the charge model [Eq. (||)] and the 
electrostatics of large LDA calculations, certain discrep- 
ancies arise in these comparisons: 

(i) The LDA calculations show that the charge is not 
a single- valued function when plotted versus the number 
of unlike nearest neighbors (Fig. Q). Although there is 
not much width to the distribution for fee alloys, there 
is a significant width for bec alloys. Also demonstrated 
by Fig. [I] is that charges in the model of Eq. (ph are 
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quantized since the number of unlike nearest neighbors 
must be an integer. The LDA calculations (particularly 
for bcc) show no such quantization. 

(ii) The slope of the q — V relation {p//R\) for bcc alloys 
(Table [|) is significantly larger in magnitude in the model 
(-0.163) than in the LDA calculations (-0.115). 

(iii) There are significant fluctuations about a linear 
q — V relation obtained by the charge model; however, 
the LDA calculation show a nearly perfect linear relation 
with no fluctuations (Fig. |J). The fluctuations of the 
charge model are especially pronounced for bcc alloys. 



B. Generalizing the Model 

The charge model of Eq. (||) is based on the obvi- 
ous chemical fact that atomic charge results from charge 
transfer, and that the latter depends on the identity of 
the neighbors, since charge transfer does not occur be- 
tween chemically equivalent sites. Thus, qi should de=. 
pend on the local environment of site i. Magri et alx3 
took first neighbors to be the "leading order" contribu- 
tion to the local environment, and for the case that Magri 
et al. treated - fee alloys - we have seen that the charge 
model provides an adequate description of electrostatics. 
However, alloys based on different lattice types can have 
different structural environments, in terms of coordina- 
tion numbers and neighbor distances: In the fee lattice 
there is a significant "gap" between the distance of the 
first coordination and that of the second. In bcc, how- 
ever, the "gap" is after the second shell. This suggests 
that one generalization of the charge model which would 
affect bcc and fee alloys differently is to allow the charges 
in the model to be dependent on more distant neighbor 
shells. Thus, instead of requiring the charges to be a 
function of the number of unlike nearest neighbors, we 
define a generalized charge model in which charges are a 
function of the number of unlike nearest neighbors on the 
first several shells s of neighbors: 



Charge vs. Generalized Number of Neighbors 



^ A s ^ [Si - Si+jJ, 



(34) 



0.20 



0.10 



0.00 



-0.10 



S -0.20 



0.20 




LDA Supercell - 256 Atoms 
-Generalized Model - 2NN 



10 



12 



0.10 



0.00 



-0.10 



-0.20 




LDA Supercell - 432 Atoms 
-Generalized Model - 3NN 



2 4 6 8 10 12 
Generalized Number of Neighbors, N 

FIG. 4. Charge versus number of generalized neighbors N 
Eq. (|3^) . Shown are the predictions of the generalized charge 
model of Eq. @ using the values of A s fit to large-unit-cell 
LDA calculations given in Eq. (^6|). Also shown are the 
charges of the LDA large-unit-cell calculations of Ref. 



intents and purposes, for s > 2 in fee and s > 3 in bcc. 
Within these ranges, we found 



Af c = 0.00745, AfVAr = 0.214, A|£ 2 /Ai' 
A*> cc = 0.00786, \ ] > cc /\ ] > cc = 0.660, 

A |cc/ A bcc = 0645) a^va 1 ^ ~ 0. 



\ fee 





(36) 
(37) 



For this generalized charge model, the charge on a site 
i is linearly proportional to the generalized number of 
neighbors, N: 



^ 1 Ai 



(35) 



(s) 

where is the number of unlike neighbors in the sth 
shell for the atom at site i. In this new "generalized" 



charge model of Eq. (34), there are S parameters, where 



S is the number of shells included. 

To determine the parameters of the generalized charge 
model, we have fit (via a least-squares procedure) the 
charges of the LDA supercell calculations to Eq. (|34| ) 
including five shells. The parameters X s are zero, for all 



C. Testing the Generalized Model 

We show in Figs. ||, |4| and |B| results for generalized 
fee (bcc) charge models including the first two (three) 
neighbors shells with these values of A s . 

Figure || shows the charge-charge correlations of ran- 
dom £=1/2 alloys predicted by the simple charge model 
of Eq. (||), and the generalized charge model of Eq. (|34|). 
Although there are currently no LDA results with which 
to compare (due to the size of the current LDA super- 
cells) , we note that the generalized charge model changes 
the sign of the second neighbor correlation in both fee and 
bcc relative to the simple model. It would be interesting 
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FIG. 5. Charge versus Coulomb shift as predicted by the generalized charge model of Eq. (E 
i?i given in Eqs. (||) and @). 



using the values of A s and 



to compare these correlations with those of LDA when 
larger supercell calculations become available. 

The charge qi versus the generalized number of neigh- 
bors N is shown in Fig. [| for LDA and for the model 
of Eq. ([54J). For fee alloys, the corrections induced by 
generalizing the charge model are small since the original 
model of Eq. (||) is already very good. The predictions of 
the generalized charge model fit the LDA supercell data 
extremely well even for bee alloys, where the nearest- 
neighbor model of Eq. (||) was lacking. 

Figure || shows the relation between charge qi and 
Coulomb shift Vi of the generalized charge model, com- 
paring the results with LDA. Like LDA, the generalized 
charge model predicts a linear relation between qi and 
Vi with almost no fluctuations. Furthermore, the slope 
of these linear relations are in excellent agreement with 
the LDA supercell data (Table ||), provided that cutoffs 
for fee and bee are at second and third neighbor shells, 
respectively. Thus, the generalized charge model of Eq. 
(|34) rectifies all of the discrepancies noted above (Sec- 
tion |v]A) between model and LDA calculations. [The fee 
model for nearest-neighbors only is already accurate with 
respect to LDA calculations (Figs. |l| and g|) , thus general- 
izing the fee charge model to first- and second-neighbors 



does not produce a large effect.] In Fig. g, we show the 
values of the parameters A s versus distance of the shell 
s. One can see that the parameters are reasonably well 
fit by an exponential function 

A1-R1 



A., 



(38) 



with a decay length of i?o = 0.34i?i. This suggests that in 
an alloy the net charge on each site is screened effectively 
in a very short range. 

Since the generalized charge model predicts a linear 
q — V relation in disordered alloys, with almost no fluc- 
tuations, one can also obtain a generalized model of the 
Coulomb shifts in an alloy 



(39) 



fc s =i 



Thus, the Coulomb shifts, like the charges, depend only 
on the occupation of the first few neighboring shells. 



D. Extracting Values of A from LDA: 
Supercell-Size-Dependence 
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FIG. 6. Charge transfer parameters \ a of generalized 
model of Eq. (BJ) as a function of distance. Also shown 
is a fit to the parameters to the exponential function of Eq. 
@. The fitted value is R = 0.34R X . 
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FIG. 7. Charge versus number of generalized neighbors N 
Eq. (^) . Shown are the predictions of the generalized charge 
model of Eq. (^) using the values of X s fit to small-unit-ccll 
data given in Eq. (^). Also shown are the charges of the 
LAPW small-unit-cell calculations of the present work. 



We have demonstrated the validity of the generalized 
model of point charges and shown how the parameters of 
the model A s may be extracted from large-unit-cell LDA 
calculations. However, the models of point charges (both 
the simple and generalized models) assume that the phys- 
ical mechanism underlying excess charge on a site is the 
same for ordered and random alloys. This suggests that 
the values of A s could be obtained from small-unit-cell 
calculations. For computational simplicity, one should 
know whether it is equally valid to extract values of A s 
from ordered or random alloys, and whether one can use 
LDA calculations of small cells (~2-4 atoms) to extract 
the values of A s . To this end, we have complementetLthc 
large unit-cell LDA calculations of Faulkner et alx3 on 
random bcc Cu-Zn alloys by performing calculations of 
several ordered small-unit-cell bcc-based Cu-Zn ordered 
compounds. We use the linearized augmented plane wave 



(LAPW) methodB The ordered structures considered 
are all bcc superlattices: CuiZni (001), Cu 2 Zn 2 (111), 
Cu 2 Zn 2 (001), Cu 2 Zm (001), Cu 3 Zm (111), and CuiZni 
(101). All of these compounds have 2-4 atoms/cell and 
the first five are commonly referred to by their Struc- 
turbereicht designations: B2, B32, BU, ClU, and DO3, 
respectively. In the LAPW calculations, we have used 
the exchange correlation of WignerE3 The muffin-tin 
radii are chosen to be equal (2.2 a.u.) for both Cu and 
Zn. Brillouin-zone integrations are performed using the 
equivalent k-point sampling methodEd. Since the charge 
model is appropriate only for charges in unrelaxed geome- 
tries at fixed volume, all computations were done in ideal 
geometries at a single volume (a=5.56 a.u.), even though 
several compositions are considered. The excess charges 
were extracted from the LAPW calculations by integrat- 
ing the charge density inside the muffin-tin spheres and 
dividing the interstitial charge evenly between the atoms 
in the unit cell. 

The LAPW charges for the six small-unit-cell com- 
pounds calculated were fit to a form of the generalized 
charge model of Eq. (p3) with lst-3rd neighbor shells. 
The parameters of the generalized model fit to these small 
unit-cell calculations, 

A^ cc = 0.00680, A^ CC /A5 CC = 0.609, Aj^/A^ = 0.131, 

(40) 

agrea-well with those fit to large unit-cell data [Eq. 
@]£3 The parameters of Eq. @ fit to small-unit-ccll 
LDA calculations lead to a q — V relation which is linear, 
with no fluctuations, and has a slope of 7/i?i=0.112, 
compared with 7/i?i=0.119 for the parameters of the 
generalized model fit to large-unit-cell LDA data. The 
charges extracted from small-unit-cell LDA calculations 
are shown in Fig. ^ as a function of generalized number 
of neighbors TV [using the values of X s in Eq. (puj) ] . These 
calculations demonstrate that the parameters of the gen- 
eralized model may be determined from calculations of 
several small-unit-cell ordered compounds in unrelaxed 
geometries at fixed volume. If one wishes to assess the 
explicit volume-dependence of the parameters, one only 
needs to repeat these types of calculations at a few differ- 
ent volumes. We have performed such volume-dependent 
calculations at a=5.36 and a— 5.75 (in addition to the 
a=5.56 calculations described above), and find that the 
values of A only have a slight volume-dependence in this 
range: The value of Ai at a=5.36 is about 6% larger 
in magnitude than Ai at a=5.75. Also, in this volume 
range the ratios A2/A1 and A3/A1 vary by less than their 
uncertainty due to the fit. 



VI. CONCLUSIONS 

Recent0Ei large scale (256-432 atom) local density ap- 
proximation (LDA) supercell calculations of Cui-^Zna, 
random alloys allow us to examine the adequacy of sim- 
ple models describing the dependence of point charges in 



12 



disordered alloys on the atomic environment. We find 
that a model in which the excess charge qi on an atom in 
an ordered or random alloy depends linearly on the num- 
ber iV,- of unlike neighbors in its first coordination shell 
correctly describes the trends in charge versus number of 
unlike nearest neighbors, (particularly for fee alloys), the 
magnitudes of Coulomb energies in random Cui-^Zn^, al- 
loys, and the relationships between constant- occupation- 
averaged charges and Coulomb shifts (V^) in the ran- 
dom alloy. However, for bee alloys the fluctuations pre- 
dicted by the model in the qi vs Vi relation exceed those 
found in the LDA supercell calculations. Although we 
found that the fluctuations present in the model have a 
vanishing contribution to the electrostatic energy, gener- 
alizing the bee (fee) model to include a dependence of the 
charge on the atoms in the first three ( two ) shells (rather 
than the first shell only) removes the fluctuations from 
the model, in complete agreement with the LDA data. 

Other possible generalizations of the charge model in- 
clude: (i) non-linearities in the charge as a function of 
number of neighbors and (ii) charges which depend not 
only on the number of nearest neighbors, but also on 
the particular arrangement of the neighbors. This type 
of dependence would lead to not only pair correlations 
among charges, but also multibody correlations. Cur- 
rently, there are no indications that these types of gen- 
eralizations are warranted. 

We thank Drs. Y. Wang andX.. M. Stocks for pro- 
viding us their LDA data of RefO. We also thank Dr. 
Z.-W. Lu for insightful discussions. This work was sup- 
ported by the Office of Energy Research (OER) [Division 
of Materials Science of the Office of Basic Energy Sciences 
(BES)], U. S. Department of Energy, under contract No. 
DE-AC36-83CH10093. 

APPENDIX A: Analytic Derivation of the q { -Vi 
Relation Within the Charge Model of Eq. (g). 

Here we derive the q — V relation predicted by the 
charge model [Eq. (^)], averaging over any fluctuations. 
Consider a random A\^ X B X alloy at x = 1/2 with nearest 
neighbor coordination Z\ and an A atom at a central 
site, denoted by A(0). (There is no loss of generality 
by choosing this atom to be A.) The charge on A(0) 
has the distribution q^j — — 2MA (M = 0, Z\) with the 
probability 



1 Zl 

(Em) = 5 £ 



puquV m 



(43) 



M=0 



where V m is the Coulomb shift on the central site, av- 
eraged over all configurations where there are M B(\) 
atoms on the nearest neighbor shell. Thus, we need to 
determine V m as a function of qu = — 2MA where 



(44) 



In order to compute Vm, we first need to compute 
q m (M). 

First Shell: For the nearest-neighbor shell, m=l, 



Ql (M) = Mqf + (Zx - M)qf 



(45) 



For the Z\ nearest neighbors of an atom in this first shell, 
one is A(Q), K\ are also nearest neighbors of ^4(0), and 
Z — Z\— K\ — 1 are remaining. For each A(l), the prob- 
ability that it has n B neighbors (i.e., with charge — 2nA) 
I of them come from atoms which are not neighbors of 
A(0) is 



Pn,l 



(M\(Z\ — 1— M\ 
1 \K ){ K x -K ) 



I J 2^1-ifi-i 



(V) 



(46) 



where K = n — I and the following inequalities must be 
satisfied 

< n < Z x - 1 
< I < n ; I < Z 
0<n~l<Ki 
n-l <M 

Kx - (n - 1) < Z x - 1 - M (47) 

Similarly, for each B(l), the probability that it has n A 
neighbors (i.e., with charge 2nA), I of them which are not 
neighbors of A(0) is: 



B(l) 

Pn,l 



Z 



2Z 1 -K 1 -X 



Zx - My M - 1 
K )\Kx-K 

(V) 



(48) 



where K = n — I — 1 and the following inequalities must 
be satisfied 



Pm 



1 

2^ 



Zx 
M 



Therefore, the energy of the random alloy is 



1 Zl 

Af=0 



PMqM 



JXrr 



q m (M) 



(41) 



(42) 



where q m (M) is the sum of charge on the mth shell sur- 
rounding A(0) under the constraint that there are M 
B(l) atoms on the nearest neighbor shell. R m is the dis- 
tance of the mth shell atom from A(Q). (Em) can also 
be written as 



l<n< Zx 
0<Z<n-l; I <Z 
0<n-l-l <K X 
n-l -I < Zx - M 
K x -{n-l-l)< M - 1 



(49) 



Combining Eqs. ( [45|fl9| ), we have 



qx(M) = (Zx-M) -2^Y,Pn 

n=0 1=0 
n=l 1=0 



A(l) 
I 



(50) 
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where p^p and p^P are subject to the constraints ( f47| ) 
and 

More Distant Neighbor Shells: For m > 1, 



g m (M) = -^fe*+g£] (51) 

Atoms on the mth shell have Z\ nearest neighbors, K m 
of them are also nearest neighbors of A(Q). Therefore, 




(52) 



where K = n — I and the following inequalities must be 
satisfied 

< n < Zt 

< I < n ; I < Zi - K m 

0<n-l<K m 

n-l <M 

K m -(n-l)<Z 1 -M (53) 

and 

( Zi-M\( M \ 
B(m) _ ( Z Y -K m \ 1 V K )\K rn -K )_ 




(54) 



where K = n — /, subject to the following constraints 
< n < Zi 

0<l <n; I <Zi- K m 
< n - I < K m 
n —1<Z\ — M 

K m -{n-l)<M (55) 

Combining Eqs. (^l||55|), we have 

1 _ Zl . " 
q m (M) = \z m £ 2n*5>3 m) - p5 m) ] (56) 

n=Q 1=0 

where and ;™' are subject to the constraints ( |5^ ) 

and (|55|). Note that q m (M) = for any shell which does 
not share nearest neighbors with A(0) (i.e., K m = 0). 
Using the above derived values of q\ (M) and q m (M) in 
Eq. (f44|), we may determine as a function of qu — 
—2MX, and as a function of M this relation is precisely 
linear, with no fluctuations: 

9m °c V M - (57) 



APPENDIX B: Structural Information for 
SQS-16 



The ideal (unrelaxed) fcc-based SQS-16 structure 
(AgBs) has orthorhombic symmetry and primitive lat- 
tice vectors 

a=(ii,0)a; b=(l,-l,2)o; c = (l,-l,-2)a 

(58) 

The 16 atomic positions, in Cartesian coordinates, are 



A : 


(0,0, 0)a 


B : 


A In 

v 2 2' 




(- - Da 

( 2' 2' 1)a 


FS ■ 


(1 -i -)« 


A : 


( 1, — , — a 
v ' 2' 2 7 


A : 


(1,-1,-1)0 


I J . 


( 2' L) 2 )a 




f 3 3 0)a 


B : 


(- - -Da 

4' 2' ' 


A : 


v ' 2' 2 y 


B : 


(1,-1, 0)a 


A : 


(2'" 1 '2 )a 


A : 


,1 In 

(2'°-2) a 


B : 


( 5 ,-2,0)a 


B : 




B : 


(1,-1,1)0. 



(59) 



The SQS-16 structure matches the first seven pair corre- 
lation functions of the random x = 1/2 alloy exactly. 
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